Electrodynamics in the presence of an axion 
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Abstract 

The low energy limit of an axion field coupled to gauge fields is 
investigated through the behaviour of the gauge fields propagator in 
a local vacuum angle background. The local (singular) part of the 
effective action for the axion field is calculated at one loop level. In the 
case of a timelike, linearly growing axion field, representing a massive 
axion, we give an asymptotic expansion of the causal propagator and 
we solve nonlocally for the first coefficient. We show that, for a generic 
axionic background, short distance propagation of the gauge field is 
well defined. 
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1 Introduction 



One of the most interesting aspects of instanton effects in QCD is related to 
the strong CP problem and why the vacuum angle is extremely small fl[]. 

Various solution have been proposed in the literature for this problem. 
In one scheme one can set 9 = at tree level and then assume that higher 
order effects in the magnitude of 9 are negligible 0, while a second option is 
to introduce a massless quark so that 9 becomes unobservable. This second 
possibility is ruled out phenomenologically. 

A third alternative is to introduce a dynamical variable 9(x) in such a 
way that an effective interaction 

L int = c9(x)FF (1) 

(c being a constant) 

dynamically sets the vacuum angle 9qcd to be zero ||. 

This third approach, generally known as the Peccei-Quinn solution, solves 
the CP problem and has cosmological implications. In this scenario 9 is 
therefore made local, where nonperturbative QCD effects produce a mass for 
the axion and fix the background value of 9qcd to be zero. 

The interaction of an axion field to electromagnetic fields was initially 
investigated several years ago, after it was shown |4[] that in a 9 vacuum a 
t'Hooft-Polyakov magnetic monopole acquires an electric charge proportional 
to 9 and to the magnetic charge g of the monopole, 

q e = C9g (2) 

with C a constant. In particular the dynamics of monopoles traversing ax- 
ionic domain walls has since been elucidated 0,0. 

In its full generality, the interaction of the axion with the electromagnetic 
field is assumed to be given, at low energies, by the action density 

L = --F 2 + -3^9 - ^-[1 - cos(Na/v)] + aFF (3) 

where F denotes the dual of the field strength and 

e 2 /. T e NB{x)\ 
a =32^K{ 6qCd+ ^v) (4) 
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where 9(x) is the axion field, m is the axion mass, N is the number of axion 
vacua and v is the vacuum expectation value that breaks the Peccei-Quinn 
symmetry. 

9° w is the electroweak angle at the grand unification mass scale (GUM) 
and 9 qc d, the constant is taken at GUM. 
Tq = 2n for QCD, is the period of 9. 

We have assumed that there is grand unification of the strong and elec- 
troweak interactions and we have used 

g s = e 2 /stn 2 (9° w ) 

at GUM, where g s is the QCD coupling constant at the same scale. 

In this letter we investigate the propagation of high frequency electromag- 
netic modes in the presence of an axion. We use the light cone expansion of 
the propagator for the electromagnetic field, and we calculate the local part 
of the one-loop effective action of the axion field. In the case of a slowly 
varying 9 field, representing a massive axion, we show that the coefficient 
of the leading singularity of the asymptotic expansion can be explicitely de- 
termined in full nonlocal form. This last result allows us to define a causal 
gauge field propagator in the presence of a massive axion, a result which 
cannot be obtained from the dispersion relation because it contains spurious 
poles. As noticed previously in the literature 0, this particular background 
configuration of the axion field (the massive axion) makes the propagation 
of the low frequency em modes unstable and a tachionic pole arises. The 
investigation of ref.|§ was carried out in the particular case of a coupled 
Maxwell-Chern-Simons theory, a system which shares the same properties of 
the massive axion model. 

We show that, for any background axion field, and in particular for a 
massive axion model, the em propagator behaves well at short distances. 
Our results are based on a direct application of the Hadamard theorem [|9] 
for equations of propagation which are diagonal in the highest derivatives. 
The short distance behaviour of the model is in fact controlled by the highest 
derivatives of the gaussian operator (which is diagonal), while the instability 
is generated by the presence of lower derivatives (coupled to an e tensor) in 
the equations of motion. 



3 



The partition function of the model can be written in the simplified form 
Z = j [de}[dA tl }e i I c+je+J ^ d4x (5) 

C = -^F 2 + l -d^ l Q + -Ax)FF (6) 

F denotes the electromagnetic field, F = \eF denotes its dual. 

We will work in the Lorentz gauge. The evaluation of the effective action 
for the #-field, functionally integrating out the gauge fields is given by 

C eff = 1 -d,9d^+ l -trHP) (7) 

where 

P = g^D - cd p 6e^d a (8) 

is the relevant operator induced by the gaussian approximation. 
The equations of motion for the gauge fields are: 

d,F»» = cd^e^d a A p (9) 

where 

F^ = d^A v - cU„ (10) 
In the Lorentz gauge we find 

uA v -k vaP d a AP = Q (11) 

where we have defined 

kvafi = cd^ aP (12) 

The Hadamard expansion || for the propagator of the gauge fields in the 
axionic background can be set directly by introducing the Green's function A 
for eq.(ll), since the gaussian operator is diagonal in the highest derivatives 

Uk vp - k ua(3 d a & p = 5 vp 8\z) (13) 

where z = x — y. The support of the A distribution is therefore inside the 
light-cone and causality is respected. 
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The standard ansatz is therefore 
A up (x,y) = Gf(x,y)D F 



oo Jl 



i 7 l °° 7 A 1 

where 

Df = 4*H(z*-iO) (15) 

is the Feynman free propagator in coordinate space while fi is a mass param- 
eter introduced in order to keep dimensionless the argument of the logarithm 
in the expansion. 

The recursion relations for the coefficients of the expansion are easily 
obtained. 

By equating to zero the independent singularities of eq. (13) applied to 
the expansion (14) we get a leading singularity equation for z~ 4 

2z^G^(x,y)-z a K aP Gf p (x,y) = (16) 
a, D F (i.e.z~ 2 ) singularity equation: 



-K aP G ( ^ p = (17) 



and a log-equation (for z 2n ln(j^)) 



z^G n v f + (n + 2)G^ + 

- y KapG^ ~ K a ^Gf p +1) = (18) 

Defining G^~^ = 0, n = —1,0, 1, we summarize all the equations in 
the form: 

z»d p G% + ^ + (n + l)G% + ^ + DGW 

- y k ^ G t 1] ~ kuapd a G p n ^ = (19) 
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At this level the axion is a background field and we see, for instance from 
eq. (16), that it affects the propagation of the leading singularity in a non 
trivial way (through k va p). 

For a constant 0-angle the equation for the propagator of the leading 
singularity is trivially given by 



tfd^l = (20) 
valid on the light cone surface, with the initial condition 

GM(x,x)=5 vp (21) 

The solution is trivially given as the Kronecker delta, and the strength 
therefore is diagonal over the entire characteristic surface. In this simplified 
case the equations of motion for the gauge fields are given by 

□ = 5 up 5 4 (z) (22) 



where as usual, 



1 = ^vp = Ovp (23) 

At this point we focus our attention on eq (16), which describes the 
leading behaviour of the dynamics of the electromagnetic field in the presence 
of a local vacuum angle. 

We will preliminarly show that, in the case of a linearly growing axion 
field, in the timelike case, some components of the leading singularity tensor 
are constant on the light cone surface. The reason appears to be quite 
simple and is valid both in the abelian and in the nonabelian case. In fact, by 
contraction of both terms of (16) with d u 6 and using a symmetry argument 
(the antisymmetry of k defined in eq (12) ) we get the equation 

d v 9z' i d li G ( ® vp (x(s),0) = (24) 

If we introduce the parametrization 

x »(s) = sx»,0< s < 1; y = (25) 
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for a straight line inside the light cone surface we can rewrite eq (24) in 
the form 

d v e(x)^°\ p (x(s),0) = (26) 
Repeating the same procedure a second time one gets the equation 

x^G^(x(s),0) = (27) 

It is simple to show that if the variation of the 9 field is linear in a timelike 
direction, a frame can be found in which it has only a time variation, say 

d„e = (ao, 0,0,0) (28) 

and in particular we get for the timelike components 

1«<? - (29, 

From the initial condition eq (21) therefore we get 

G (0 \ t = G (0 \o = l (30) 

As we emphasized previously, this result remains true even in the more 
realistic non abelian case, for a linearly growing (timelike) axion field, and 
is a simple consequence of the antisymmetry involved in the coupling. A 
simple and physical way to look at this result is to view it as an anisotropy 
effect induced on the propagation of the gauge fields by the pseudo tensor 
coupling. 

We will show now that the other components of the leading singularity 
tensor can also be determined by following an approach analogous to the one 
we used above, in the case in which condition (28) is enforced. 

The equation for the leading singularity (eq. 16) can be cast in the form 



_G«\.(x( s ),0) - ^^e ^G(°)^(x( S ),0) = k, I = 1,2,3 (31) 

For a given timelike vector x^, the previous equation has two first inte- 
grals. This can be easily seen through an elementary analogy with the motion 
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of a classical charged particle under a Lorentz force in a constant magnetic 
field. We use this analogy to solve it. Define 

Gf(x(s),0) = v f , u k = ^x k (32) 

where Cj is constant at fixed x^ and rewrite it as 

jfj = & A Vj (33) 
and two first integrals of motion along the s-line are 

Cj.v j = ^j*G^(x(s),Q) = d s (34) 

v 2 j = E G^(x(s), 0)G\f(x(s), 0) = l 3 (35) 

i 

Iterating eq.(33) we get 

By usingff the initial condition eq (21) one easily gets 

dj = ^x j lj = l, J = 1,2,3 (37) 

and finally the solution 

Gif{x{s), 0) = [dij - ^-^J cosi^-rs) + e^ysm(yrs) + (38) 

This expression describes the behavior of the gauge field correlator in its non 
local form, asympotically around the light-cone. 

Our results can be summarized in the following expressions for the Feyn- 
man propagators of electromagnetic fields in the linearly growing timelike 9 
vacuum 
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1 

+ /og. terms 



47r 2 i(x 2 — iO) 

(39) 



< 4,(3)4(0) = o) + Gfffo 0)%(a; 2 - iO) (40) 

modulo additional logarithmic corrections. This model provides an addi- 
tional example of a solution for the Green function in an external field in full 
nonlocal form. The outline of this approach, as is well known, dates back to 
Schwinger [|11|] and it dealt with the Dirac propagator in an abelian gauge 



field. Causality problems have been investigated by Velo and Zwanziger |10] 



in this framework successfully, and calculations of anomalies are also possible. 
The massive axion model suffers, however, of an instability at low frequencies, 
in other words the retarded em propagator developes a tachionic behaviour 
in its lower modes. This aspect has been considered in very detail in ref. 
[§] in which the 3+1 dimensional Maxwell-Chern-Simons (MCS) system is 
discussed. The equations of motion for the propagation of electromagnetic 
fields in the massive axion background are equivalent, in fact, to the MCS 
model. An expression of the retarded propagator was also given in it and it 
was shown that em frequencies smaller than do, the time component of the 
derivative of the axion field, becomes tachionic |J. We refer to that paper 
for further details. In the MCS case the instability is an inevitable effect, 
since no background field is involved from the beginning, while in the massive 
axion model is still unclear whether by taking into account in a better way 
the axion dynamics the instability can be removed. 

In its current formulation, in the case of either a spacelike or timelike 
variation of the local vacuum, the propagation of the gauge fields appears 
to be consistent in the ultraviolet regime and, therefore, the singular part of 
the effective action for the axion field, obtained by integrating out the gauge 
fields in the partition function, can be computed straightforwardly. 

In the following we are going to show that the structure of the singularities 
of the propagator of the gauge field in an arbitrary 9 background has a very 
simple feature. 
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2 The effective action 



In the general case of a variable vacuum angle we can easily calculate the 
singular (local) part of the effective action for the 9 field integrating out the 
gauge field. Let's define [G^J] to be the coincidence limit of the coefficients 
introduced before. 

[g; u ] = g$(x,x) (4i) 

The same coefficients also appear in the divergent part of the one loop effec- 
tive action 

rW(0) = ^rm(P) (42) 
Let's take a variation of (42) with respect to the background field 9 

5T^[9] = UrdPP- 1 = % -< 5P, u A^(x,y)5 4 (z) > (43) 

where the angular bracket denotes integration over x,y. 
We find after some algebra 

[5P vP k^] = -c< d,89e^ av d a A vp {x, y)5\z) > (44) 

Ultraviolet divergences in the effective action come from the singular be- 
haviour of the propagator at short distances. We need to regularize expres- 
sions such as [ln(^ 2 — iO)], [D F ] and [d^D F ]. This can be done in many ways. 
The simplest one consists in introducing a short distance cutoff A. All the 
singularities are regulated at the intermediate stage, while we first compute 
the coincidence limits of the various coefficients. The cutoff is removed by 
taking the limit A — > oo at the end. 

A simple calculation shows that the relevant contribution to (44) is given 

by 

[d a Kvp] = -^M* 2 " iO)][9«Gg)] + [d a G<$][D F ] + [G^ up ][d a D F ] (45) 

For notational simplicity we omitted to write down in eq. (45) explicitely 
the cutoff dependence of the singular coincidence limits. 

It can be seen immediately that, in this regularization, the contribution 
from the [dD F ] singularity is zero. 
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By differentiating the singularities equations we find 

[d m G^] = \k ump (46) 
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g kvm^kflnp ~l" gkvn^kflmp (^7) 



The explicit expression for the counterterm is 

< [8P TV A"] >= J d 4 x5d,6e^ (^r[ln(^ 2 - iO)][d a G$] + [d a G<® vp ][D F \ 

(48) 

A straightforward but lenghty calculation gives 



r<1, «" = U** & n ^l^ 2 - + m&T? 

(49) 

where we've defined 

Q 2 = d^O (50) 

Introducing a pole and a logarithmic renormalization costants Z\ and Z2, we 
can express the counterterm lagrangean in the form 

VC count. = cZi q2 + c ^ ( C 2(Q2)2 _ 2d /1 0Dd^) (51) 

3 Conclusions 

We have analyzed the dynamics of a coupled axion-gauge field model in the 
low energy limit and we've investigated in detail the problem of the prop- 
agation of high frequency em modes in an axion background. In particular 
we have seen that, in the case of a timelike varying axion background, the 
leading singularity of the Hadamard expansion of the gauge fields correlator 
can be determined in non local form. By applying Hadamard's theorem on 
the structure of the singularities of diagonal hyperbolic operators, we have 
shown that the instability of the Maxwell-axion system is confined to the 
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propagation of low em frequencies, reobtaining a result first considered in 
ref. ||. Our results are valid for the propagation of em modes in any axionic 
background. We have also shown that the local part of the effective action 
of the axion can be computed in a simple form, a result which gives to the 
effective model physical consistency. 
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